Transport in finite graphene samples with random gap 
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We study the DC transport of finite graphene samples with random gap. Using Dirac fermions to 
describe the low-energy physics near the Dirac point, we employ a generalized Drude form for the 
conductivity. The latter is constant for a vanishing average gap but always decreases with increasing 
sample size for a nonzero average gap. The asymptotic conductivity of the infinite sample is either 
nonzero if the average gap is smaller than a critical value or zero otherwise. Our results are in 
qualitative agreement with recent numerical calculations of Bardarson et al., Phys. Rev. B 81, 
121414(R) (2010). 
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Transport in graphene is a fascinating field in 
terms of experiments as well as theory [l|-|4( and it 
provides an interesting test ground for the validity 
of general theoretical concepts. A problem is that 
experimental measurements and numerical simu- 
lations of transport properties are performed on 
finite size samples. On the other hand, most ana- 
lytic calculations are based on linear-response the- 
ory for infinite samples. It is known from pertur- 
bative renormalization group calculations that dis- 
order in graphene can create its own finite length 
scale that describes a crossover from small scales 
(weak disorder regime) to large scales (strong dis- 
order regime) p| |g| . The perturbative renormal- 
ization group approach works only in the regime of 
weak disorder (i.e. on small scales), such that the 
crossover behavior requires an alternative analysis. 
An example is the random-gap model that is de- 
scribed by Dirac fermions with random mass. The 
renormalization of the average gap parameter rh 
always scales to zero 0] , indicating the total sup- 
pression of the gap by fluctuations. Provided that 
the perturbative renormalization group approach 
is valid, its results imply that a random gap always 
scales to zero and graphene remains metallic for 
any strength of gap fluctuations. This, of course, 
cannot be true for a sufficiently large rh, where 
we expect an insulating behavior. The latter was 
observed experimentally in case of hydrogenated 
graphene @, in a recent numerical simulation [ll| 
and analytically for an infinite random-gap model 
Q. However, to compare the analytic theory with 
the experiment and with numerical simulations 
[TTI j , the finite-size crossover remains an open prob- 
lem. In particular, the transition between metal- 
lic and insulating behavior may be washed out 
by finite-size effects. An interesting result of the 
numerical calculation is that graphene is metal- 
lic with constant DC conductivity <jq = e 2 /irh for 
vanishing average gap m, regardless of the strength 



of the gap fluctuations, but it decreases with sam- 
ple size for any nonzero average gap rh. This rises 
the question whether or not the conductivity van- 
ishes for the infinite sample with any rh > 0. This 
question will be addressed in the following. 

We employ the linear response approach to 
study the transport properties for the Dirac Hamil- 
tonian H with random gap term: 

H = iVi(7i + iV2<72 + TO r (73 . 

Oj is a Pauli matrix and the random gap m r has 
an uncorrelated Gaussian distribution with mean 
rh and variance g. Then the conductivity can be 
evaluated from the Kubo formula. The latter is 
applicable also to finite samples, provided that the 
boundary conditions are properly taken into ac- 
count [12, HH • Then the conductivity is related to 
the average two-particle Green's function as in the 
infinite sample. The average conductivity at the 
Dirac point becomes the simple scaling form [9] 
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where w is the frequency of the external field. It 
should be noticed that the conductivity depends 
on the fluctuations of the random gap only through 
the scattering rate r/. This is a generalized Drude 
formula with the scaling function F(z). The con- 
ventional Drude formula has F(z) cx 1/z. In case 
of graphene the DC limit u> = of the conductivity 
has the scaling form [9j 
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The scattering rate r\ is real, and the conductivity 
vanishes for 77 — > 0. A vanishing scattering rate 
can be understood here as the absence of quantum 
states, which reflects the appearence of an effective 
gap. 
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The scaling form of Eq. ^ allows us to define 
a critical gap parameter as m c — \f 4r/ 2 + fh 2 
which itself depends on fh. Then the conductivity 
vanishes as one approaches m = m c . In the infinite 
sample the scattering rate 77 is determined by Q 
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where g is the variance of the fluctuating gap. 77 is 
automatically zero if there is no real solution 77 for 
Eq. © [II. 

The scaling relation in Eq. (JXJ) and a nonzero DC 
conductivity are consequences of a scale-invariant 
behavior of the system due to a massless mode 
S El ■ The latter reflects a diffusive behavior of 
the quasiparticlcs, where the diffusion coefficient 
is finite and proportional to I/77 in the model un- 
der consideration. Diffusion is also possible in a 
finite (confined) system with the restriction that 
the diffusion time is limited by the sample size. 
However, the restriction can be avoided by apply- 
ing periodic or quasiperiodic boundary conditions, 
where the diffusing particles can leave the sample 
on one side and re-enter it from another side. This 
corresponds with a physical system in an electric 
circuit as it is described by the Kubo formula. 

Our previous calculations can be extended to 
finite samples with N x N sites. The contin- 
uous wavevectors of the Dirac fermions kj are 
now discrete variables with kj = 2irrij/(N + I) 
(rij = 0,1,..., N). In terms of the conductivity this 
affects the scattering rate 77. The latter becomes 
t]n in the finite sample and is determined by the 
equation 
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For a fixed parameter 77^+ to 2 /4 the sum decreases 
with increasing N. This means that for given av- 
erage gap 777. and variance g the scattering rate tjn 
decreases with increasing N. Consequently, the 
DC conductivity cto of Eq. ([2]) also decreases with 
increasing size N. If <tq > for a given finite N 
the conductivity can vanish for a sufficiently large 
sample size. Whether or not it vanishes for the in- 
finite sample depends on the asymptotic value 77, 
which is determined by Eq. ©. Thus the transi- 
tion from metallic to insulating behavior depends 
on N, where the ability to conduct decreases with 
increasing N. A special case is 777 = 0, where we 
have Co = e 2 /^h for any N, since r]pf does not 
enter the expression of Co in Eq. ©. This is in 
agreement with the results of the numerical calcu- 
lation of Ref. [Till - 
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FIG. 1: Critical gap parameter m c = -J Arjj^ + fh 2 in 
units of the electronic hopping rate of finite N x N 
samples as a function of TV for average gap fh = 0.06 
and for g — 2, 4, 5, 10. For fixed N the critical gap pa- 
rameter m c increases with increasing disorder strength 
9- 



The behavior of the critical gap parameter m c 
and the conductivity <tq is presented for sample 
sizes N = 20, 100 in Fig. [T]and Fig. H respec- 
tively. These results can be compared with the 
corresponding results for the infinite sample given 
in Table fl] All calculations use the same average 
gap ?T7 = 0.06 for easy comparison. 

There are three remarkable facts for fh > 0: (i) 
The conductivity and the scattering rate always 
decay monotonically with size N. (ii) The slope of 
the nonzero conductivity decreases with increasing 
variance g. In particular, the conductivity is al- 
most constant for g = 10 (cf. Fig. [2]). (iii) At fixed 
N the conductivity and the scattering rate increase 
always with increasing disorder g. These results for 
the conductivity also agree with the findings of the 
recent numerical simulation of finite samples with 
random gap (Tlj . although a quantitative compar- 
ison is not possible due to different types of gap 
randomness: Instead of the single parameter g of 
uncorrelated randomness in Eq. ((4]), the numerical 
simulation uses correlated randomness with corre- 
lator g r - r '- The identification of g with the sum 
Kq — ^2 r g r may overestimate the strength of ran- 
domness in case of correlated disorder [15| . More- 
over, our analytic result allows us to extrapolate 
to the infinite sample. This reveals an insulating 
(gapped) behavior for fh > m c but also a metallic 
behavior for fh < m c , where m c is determined for 
the infinite sample (cf. Table HJ. 

The size dependence of the conductivity is sim- 
ilar if we send the width of the sample to infinity 
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FIG. 2: Minimal conductivity ao of finite N x N sam- 
ples as a function of N for average gap m = 0.06 and 
g — 2,4, 5, 10. For fixed N the conductivity increases 
with increasing disorder strength g. 
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m c 


a (e 2 /nh) 


2 


0.06 





4 


0.06 





5 


0.09 


0.78 


10 


1.12 


0.98 



TABLE I: Critical gap parameter and minimal conduc- 
tivity of the infinite sample for average gap fh = 0.06. 
(For this calculation the area of integration is approx- 
imated in Eq. © by k 2 < I.) 



and study a variable sample length L. Then one of 
the sums in Eq. (U) becomes an integral that gives 

1 A arctan(27r/ ^ff + 4n 2 n 2 / (L + l) 2 ) _ 1 
L + l ^ 27iyij 2 +4^2/(^ + 1)2 ~ g ' 

The g and L (or N) dependence for the N x N 
square sample and the infinite strip is plotted in 
Fig. 1 

Our observations in (i) - (iii) provide a quali- 
tative interpretation of the transport properties in 
finite graphene samples in the presence of a ran- 
dom gap. First of all, the expression of Co in Eq. 
(P) implies that the conductivity depends on the 
random gap only through the average gap fh and 
the scattering rate r\. rj itself is determined by Eq. 
Q (infinite sample) or by Eq. ^ (finite sample), 
respectively, as a function of fh and of the variance 
of the gap fluctuations g. There is no effective gap 
as long as the scattering rate is positive, since the 
conductivity is nonzero. The fact that at fixed fh 
the scattering rate increases with increasing g is 
indicative of an increasing density of states near 
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FIG. 3: Conductivity as a function of disorder strength 
g for a fixed size N x N with N = 40, 50, 60, 80, 100. 
Larger values of iV require larger values of g for the 
same conductivity. 
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FIG. 4: Conductivity as a function of disorder 
strength g for infinite width and finite lengths L = 
40, 50, 60, 80, 100. Larger values of L require larger val- 
ues of g for the same conductivity. 



the Dirac point due to disorder, since gap fluctua- 
tions allow the appearence of additional states. It 
is less obvious whether these states are localized 
or extended, although the increasing conductivity 
with increasing g supports the latter. The con- 
stant conductivity for fh — can be understood 
as a balanced interplay of an increasing density of 
states near the Dirac point and a decreasing diffu- 
sion coefficient due to gap fluctuations. Our scal- 
ing relation in Eq. @ yields a linear decay of the 
conductivity when m c is approached. This is in 
agreement with a recent numerical simulation for 
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Dirac fermions with a uncorrelated random mass 
(l6j |. It remains an open question whether or not 
the linear behavior of the conductivity near m c is 
universal with respect to different types of disor- 
der or other details of the model on short scales. 
Numerical simulations for a similiar model have 
revealed that the conductivity obeys a power law 
near m c with a different exponent i/»l,4 [l7j . 

In conclusion, we have found that transport in 
finite graphene samples with random gap is en- 
hanced by the gap fluctuations. In particular, the 
DC conductivity at fixed average gap increases 



with increasing variance. The samples are metallic 
for sufficiently small but nonzero average gap. On 
the other hand, the conductivity always decreases 
with increasing sample size. 
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